Fidelity for kicked atoms with gravity near a quantum resonance 
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Kicked atoms under a constant Stark or gravity field are investigated for exper- 
imental setups with cold and ultra cold atoms. The parametric stability of the 
quantum dynamics is studied using the fidelity. In the case of a quantum resonance, 
it is shown that the behavior of the fidelity depends on arithmetic properties of the 
gravity parameter. Close to a quantum resonance, the long time asymptotics of the 
fidelity is studied by means of a pseudo-classical approximation first introduced by 
Fishman et al. [J. Stat. Phys. 110, 911 (2003)]. The long-time decay of fidelity 
arises from the tunneling out of pseudo-classical stable islands, and a simple ansatz 
is proposed which satisfactorily reproduces the main features observed in numerical 
simulations. 
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I. INTRODUCTION 

The stability of quantum evolution against parametric changes of the quantum Hamil- 
tonian is a subject of wide theoretical and experimental interest. A widely used concept 
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is the fidelity introduced by Peres and the closely related Loschmidt echo [2J, y|, which 
is built as the interference pattern between states which are obtained by propagating the 
same initial state under Hamiltonians, say Hq and H, which are slight perturbations of each 
other. A standard definition of the fidelity is: 
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F(t) 



i e iH t/h e -iHt/h\ 



(1) 



The behavior of fidelity in time is known to display some universal properties that reflect the 
underlying classical dynamics {2, 3]. Such properties have been mostly explored for the case 
of systems, which are chaotic in the classical limit. In this paper we will study the fidelity 
in the mixed phase space regime. For this the system under consideration is the quantum 
kicked rotor 4|-[6]. 

The main motivation for our analysis here is twofold. Firstly, there has been a growing 
interest over the last decade in the dynamics of the quantum kicked rotor (and its variants) 



at and close to the so-called quantum resonances [7j, theoretically (see, e.g, |s-16|) as well 
as experimentally (see, e.g., |17H23|). Secondly, only recently concepts have been developed 
to actually access the fidelity in setups based on cold or ultracold quantum gases. The used 



techniques range from interferometric methods in either internal atomic states [24| or in the 
center-of-mass motion of the atoms 23[ to the time reversal of the dynamics by exploiting 
the properties of the quantum resonant motion 2jJ, 125 ]. 

In this paper we studythe quantum kicked rotor under the additional influence of a 
Stark or gravity field In Sect. [TT] the Hamiltonian of one kicked atom is quickly 

reminded and the fidelity, which is the main quantity studied here, is precisely defined 
for our system. We then report on a subtle dependence of the fidelity on the arithmetic 
properties of the relevant parameters (Sect. ITH|) . a result which may be interesting for future 



precise measurements of fundamental constants (see the discussions in refs. 
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23, 



26|). 



Sect. [IV] is devoted to the dynamics close to quantum resonance. Based on the pseudo or 
e— classical formalism developed by Fishman et al. |8|, we explain the overall behavior of the 
fidelity using e— classical phase space densities and quantum tunneling rates from the stable 
resonance island to the surrounding chaotic sea in phase space. Some technical details are 
found in appendices. 
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II. THE KICKED ROTOR WITH GRAVITY 

We are interested in the quantum dynamics of a particle moving in a line, periodically 
kicked in time, and subject to a constant Stark or gravity field. It is described by the 
Hamiltonian in dimensionless variables (such that h = 1) jsj]: 

"2 

H(x,p,t) = ^- - - x + k cos(x) V s S(t - tr) . (2) 

tez 

The kicking period is r, the kicking strength is k, and t is a discrete time variable that 
counts the number of kicks. The parameter r\ yields the change in momentum produced 
by the constant field in one kicking period. In the accelerated frame of reference [sj], the 
potential experienced by the particle is periodic in space and so the quasi-momentum (3 is 
conserved by the evolution. With the chosen units, takes all values between and 1. Using 
Bloch theory, the particle dynamics can then be identified with that of a family of quantum 
rotors, labelled by the values of (3. For the /3-rotor (i.e., the rotor in the family to which a 
given value (3 of the quasi-momentum is affixed) the evolution from immediately after the 
(t — l)-th kick to immediately after the t-th kick is described by the unitary propagator f|: 

= e-' lkCOs{§) e --/2(A^+r,t+r,/2)2 ^ (3) 

and the evolution operator over the first t kicks is: 

U} AV = U PAl] (t-l) ZY AM (t-2)...ZY WjJ7 (l) U PAv (0). (4) 
where J\ is the momentum operator: 

with periodic boundary conditions. The time-dependent Hamiltonian that generates the 
quantum evolution corresponding to (BJ is then: 

^ * " ' ■ h /•' <-OS( r,' ) 

tez 

We will study the fidelity that measures the stability of the evolution (TjJ with respect to 
changes of the parameter k. For a given /3-rotor this fidelity is defined by: 



U{Af, 6,P,t) = -(tf + (3 + -t) 2 + k cos(^) S (t ~ tr) . (5) 



F /3 (k 1 ,k 2 ,ri,t) = (u} tkurj i)Up M ^) , (6) 
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Moreover, having in mind experimental situations with cold atoms [9|,lll|,H8|,|23|,l24j, we will 

also consider the case when the initial state of the atomic cloud is an incoherent mixture 

of plane waves with a distribution p(/3) of the quasi-momentum. In this case the fidelity is 
I I 

given by [121 ] : 



F(k 1 ,k 2 ,r],t) 



p{fi) (Uj M>n il> 



(7) 



III. FIDELITY AT A QUANTUM RESONANCE 



In the gravity free case (77 = 0), Eq. (pi) describes the standard Kicked- Rotor (KR) 
dynamics, and so-called KR resonances 0, 13, 2^] occur whenever r is commensurate to 
2tt. Then, for special values of quasi-momentum j3, the energy of the /3-rotor asymptotically 
increases quadratically as t — > 00. In the presence of gravity, asymptotic quadratic growth 
of energy at certain values of f3 is still possible 28]. Here we study the behavior of fidelity in 
the presence of gravity and for the case of a main KR resonance, i.e., r = 2ttI (with integer 
/). Denoting ipp(t) = Up k i/)p(0), one may explicitly compute 



(6\if}p(t)) = e --MP,vX) e -ikA(e,p, v ,t)^ e _ t ( 2j3 + 1 ^ 7r/ _ nl7]t 2 J ^(0)) 

where a((3,r},t) is a global phase, and 

t-i 

A{6, f3,t) = Y^ cos(9 - (2/3 + l)vr/ r - 2nlrr]t + nlr]r 2 ) . 

r=0 

From now on we assume that the initial state is a plane wave: |^(0)) = \n ), i.e. 



(8) 



(9) 



<0|^(O)> = (9\n ) 



2n 



(10) 



The fidelity is directly obtained from (jSJ) using the method described in 12] . First write 
A(G,P,t) = Re (e w W t ) with 



t-i 



W t = W t (v,P) = J2e~ inl{2ni)r e~ i2nlrvt ~ 



iirlrjr 



(11) 



r=0 



Then from (JSJ) it follows that the wave function after the t-th kick is given in momentum 
representation by: 



n 



3 niarg(W^-) 



1 

2^ 



2tt 



e ni ^ w t\-\y»- n j no _ n (k\w t \) 



(12) 
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where J no _„(.) is the Bessel function of order n 
formula of Bessel functions, see e.g. 7.15.(31) in 



n. Using this along with the addition 



29l | . one finds: 



(13) 



where we defined the perturbation parameter Ak = k 2 — k±, and so the fidelity of a single 
j3— rotor is given by: 



Fp{k h k 2 ,r],t) = J (AA:|Wt|) 



(14) 



A. Asymptotics of the fidelity for one single rotor 

From Eqs. ( IT2l) and (Till) it is clear that the long time asymptotics of the wave-packet 
propagation, and of the fidelity as well, are determined by the behavior of |Wt| as t — > oo. 
One may write: 

W t = e m ^W(P, V ,t), (15) 
where $(77, (3, t) = — nl(2f3 + 1 + r]t)t, and 2U(?7, 0, t) is a quadratic Weyl sum: 



i7rZ(2/3+l) r ^ml-q r 2 



(16) 



r=l 



The asymptotic behavior of such sums as t — > oo is known to depend on the arithmetic 
nature of the number 77, i.e., on whether it is rational or irrational, and in the latter case 
on its Diophantine properties 30]. First of all, as the behavior of Weyl sums may be quite 
erratic, we resort to the time-averaged fidelity defined by: 

1 

(F^h,k 2 , V ,t)) T = -J2 F dh,k2,V,t), T»l, (17) 

t=o 

which has a smoother dependence on time than the original fidelity. 

The easiest case is when rj is rational: f] = p/q, with p and q mutually prime integers. In 
that case, setting r = 2jq-\-v in the sum ( ITBl with j a non negative integer and < v < 2q— 1, 
the sum may be rewritten in the form: 



= C(r),/3,t)B(r),P,t), 
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C(r),(3,t) = -~ i,:r ' ■"' 



where 

[t,2g] 

e 
3=0 

{t,2q} 

B(v,P,t) = Yl e inl ^ +1)u e inlu2p/q , (18) 

having denoted by [t,2g] the integer part of t/(2q) and {t, 2q} = t mod(2g). The factor 
B(j],f3,t) is a periodic function of t with period 2q. Explicit calculation of the sum on the 
right hand side of the 1st equation shows that C(r),(3,t) is: 

• quasi-periodic for /3 irrational, 

• periodic for rational and 2f3q non-integer, 

• linear, i.e. C(j],f3,t) = t, when 2j3q is an integer. 

Such facts have the following implications on wave packet dynamics on the one hand and on 
the behavior of the fidelity on the other hand. Two cases have to be distinguished, according 
to whether 2(3q is integer, or not. In the former case, a quantum resonance occurs. Indeed, 
using that |Wt| is equal to t times a periodic function of t, from Eq. ( fl2l) and from the 
well-known asymptotics of the Bessel functions at large argument and fixed order, see e.g. 



7.13.1 (3) in [29|: 




/ 2 / nn 7r\ 
J n {x) ~ \ — cos I x I , x — > oo , (19) 



7TX \ 2 4, 

we find that the probability in the n-th momentum eigenstate decays in time like 1/t. Hence 
the wave packet spreads linearly in time in momentum space, and the energy quadratically 
increases. Instead, if 2/3g is not an integer then Eq. ffl2|) shows that the amplitude of the 
evolving wave function in any momentum eigenstate oscillates quasi-periodically in time, so 
no unbounded spreading in momentum space occurs. 

A similar reasoning based on Eqs. (1141) and (I17p shows that, for all nonresonant values 
of /3, the time-averaged fidelity saturates to a nonzero value in the limit T — > oo (Fig. [Th). 
Instead, at resonance (2/3q =integer) from Eqs. f|T4l , ffl9|) (with n = 0) and (fTTl) . the fidelity 
is seen to asymptotically decay to zero like log(T)/T (case (3 = 0.5 in Fig. [Ih.). 



The case of irrational rj is much more difficu 
depends on the Diophantine properties of rj 



t, as the behavior of Gauss sums crucially 
3011 . Here we limit ourselves to a heuristic 
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analysis. For strongly irrational r\ one may naively picture Wx as a sort of random walk in the 
complex plane, suggesting that | | grows like \ft in some average sense. Thanks to ffl~4]) and 
(fl9l) an asymptotic decay of the average fidelity like l/\/t is then expected. This is roughly 
numerically confirmed for the case when 77 is equal to the golden ratio in Fig. [lb. However 
the actual decay displays strong fluctuations because it depends on the continued fraction 
expansion of 77, notably large partial quotients in the latter may cause fidelity to behave as 
in cases of rational 77 over significant time scales, e.g., for 77 = n = [3, 7, 15, 1, 292, 1, 1, . . .] 
in Fig. [Tb. 
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FIG. 1: (Color online) Fidelity for the kicked rotor with gravity at the quantum resonance r = 2ir 
for k\ = 0.87T and k<i = 0.77T. The initial state is a plane wave with zero momentum no = 0. 
Several choices of 77 and /3 are shown: a) r] = 0.1. Full black line: const. ft. Blue dashed line: 
P = 0.23. Green double dot-dashed line: = 0.499. Red dot-dashed line: /3 = 0.5. b) Full black 
line: const. /vt. Blue dashed line: rj = (\/5 — l)/2 (the Golden Ratio), f3 = 0.23. Green double 
dot-dashed line: rj = it, j3 = 0.23. Red dot-dashed line: 77 = n, (3 = 0.5. 

In order to mimic the experimental setups based on cold atoms (9), [ll|, 12, 18, 23, 24], we 
consider the case when the initial state of the kicked atoms is an incoherent mixture of plane 
waves with a uniform density of 0: p((3) = 1. The expression ([7]) is computed as an average 
over a large number of randomly chosen values of 0. The result does not vary significantly 
when the number of values of /3 exceeds a few thousands. We observe a sharp difference in 
the asymptotic regime depending on whether 77 is rational or not, see Fig. [2J We again show 
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the time-averaged quantity: 

1 T_1 

(F(k h fc 2 , V , t)) T = - F (hi k * V, t) • (20) 

t=o 

For rational values of r), the fidelity is observed to saturate towards a finite value. This is 
not surprising, because this is precisely the expected behavior for all values of (3 in a set 
of full measure. On the contrary, the fidelity decays like 1/T for irrational rj, see Fig. |2j 
This is roughly explained noting that, besides the decaying prefactor ~ l^tl -1 ^ 2 i n the 
asymptotics ffl9|) (with n — 0), one more mechanism of decay is introduced by ensemble 
averaging, which affects the rapidly oscillating part of the Bessel function. The observation 




FIG. 2: (Color online) Fidelity for an ensemble of kicked rotors with gravity at r = 2tt for k\ = 0.87T 
and &2 = 0.77T. The initial state is a plane wave with uq = 0. Several choices of r\ are shown. The 
integral over f3 in Eq. ([7]) is computed via a Riemannian sum with 5000 randomly chosen /3s. Full 
black line: const. /t. Blue dashed line: r/ =Golden Ratio. Green double dot-dashed line: r/ = ir. 
Red dot-dashed line: 77 = 0.1. 



of the fine details of the number theoretical properties of r/ is certainly of high interest for 
precision measurements 2lJ, |23|, |25|, |26( but, at the same time, a big challenge to experimental 
resolution in current setups. 
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IV. FIDELITY CLOSE TO A QUANTUM RESONANCE. 

A. Reminder of the e-semiclassics. Fixed points. 

When the kicking period is close to a quantum resonance, 

r = 2ttI + e, / integer , |e| < 1, (21) 



we will implement a technique of quasi-classical approximation originally described in [8j and 
therein termed "e-classical approximation". Introducing a rescaled kick strength k = \e\k, 
and a new momentum operator 

f=\e\tf=-i\e\±, (22) 
the propagator (j3J) may be rewritten as: 

^ )M (t)=exp^-^cos(0)Jexp^-j^^(/,t)J , (23) 



where: 



Hp(I,t) = ^sgn(e)/ 2 + i[ 7T l + T {/3 + r ] t+ 1 l)]. (24) 



The small parameter e plays the formal role of a Planck constant in ff23l . so, close to 
a quantum resonance, the quantum dynamics mirrors an "e-classical" dynamics, which is 
immediately inferred from (T2"3"|) and ( 12"4"|) . After changing the e-classical momentum variable 
from I to J given by: 

J = I + sgn(e) [vrZ + r/3 + rrjt + tt?/2] , (25) 

the e-classical dynamics is described by the following map that relates the variables J and 
9 from immediately after the t-th kick to immediately after the (t + l)-th one: 

Jt+i = Jt + fcsin(0 t +i) + sgn(e)r7? 
^t+i = h + sgn(e) J t mod 2vr 

If considered on the 2-torus, this map (|26|) has a fixed point at J = and 6 = 6q if: 

sin 6» = -sgn(e)^- , (27) 
k 

and this fixed point is stable if and only if: 

< k\ cos^ol < 4 and cos9 = — sgn(e)| cos^ol . (28) 



10 



Such stable fixed points give rise to stable islands immersed in a chaotic sea. From (|25|) 
it is immediately seen that, in physical momentum space, such islands travel at constant 
velocity: 



As islands trap some of the partic 



quantum accelerator modes 



17 



v = -^. (29) 
e 

e's wave packet, they give rise to experimentally observable 



20] . Because of such modes a quadratic growth of energy 



is observed over significant time scales, both in the falling frame and in the laboratory 
frame [8|. However such quadratic growth eventually comes to an end as the modes decay, 
due to tunneling out of the stable islands Il6fl. Smaller accelerating islands may also exist, 

n n 

associated with higher-order fixed points of map (|2T)|) |14j . |15j, however we will restrict 
ourselves to the above described ones. 



B. Long time asymptotics of the fidelity 

Typical numerical results illustrating the time dependence of fidelity (jH]) are shown in 
Figs. |3H6j In all those simulations the parameters were chosen in ranges where significant, 
experimentally detectable accelerator modes exist (see Appendix [A] for details). In general, a 
very short initial transient is observed (typically up to one or a few hundred kicks depending 
on parameters), marked by a very quick drop. A clear, relatively long exponential decay 
follows. This is sometimes followed by yet another stage of exponential decay, at a slower 
rate and with stronger fluctuations. This general behavior is qualitatively understood as 
follows. The initial sharp decay is due to the part of the initial wave packet that lies in the 
chaotic component of either of the two dynamics (defined by the two different kick strengths), 
and is rapidly carried away. The fidelity is thereafter dominated by the parts of the wave 
packet which are trapped inside the islands. For this reason, in order to bring this stage of 
fidelity decay into full light, we choose our initial state in the form of a Gaussian state mostly 
located inside one island. The islands which correspond to the ki and to the k^ dynamics 
are slightly displaced with respect to each other, however they travel in momentum space 
with the same velocity ( )29|) . Therefore, the mismatch between the k\ and the &2 dynamics 
is, in e-classical terms, mostly produced by (i) different structures inside the islands and (ii) 
the decay from the islands into the chaotic sea due to dynamical tunneling. Concerning (i), 
the different rotational frequencies in the two islands are expected to produce quasi-periodic 
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oscillations of the fidelity [10|, [3l| , so (ii) should be the main mechanism responsible of the 
mean fidelity decay. This leads to the following crude description. The main contribution 
to fidelity comes from the part of each factor in the scalar product in ([6]) which is trapped 
in the respective travelling island. Hence the decay of fidelity is determined by the decay 
of each part, which is in turn determined by its respective tunneling rate into the chaotic 
sea, which we will denote by Tj for i = 1, 2. Then a simple, self-explanatory ansatz for the 
asymptotic decay of fidelity is 

F(t) oc n(Ai\A 2 ) e~ rit + fi(A 2 \Ax) e ~ r ' 2t + MA n Ax) e - (ri+r2)t , (30) 

where /x is the classical invariant measure of phase space sets, Ai is the island around the fixed 
point associated with fcj. This ansatz was found to satisfactorily reproduce the actual decay 
of fidelity in our numerical checks. In our simulations, the phase-space areas appearing 
in ( 130]) were estimated as described in appendix [B] Quantum decay rates were found as 
follows: for both the k\ and the h% dynamics we numerically computed the probability in a 
momentum range centered on the accelerator mode. This measures the amount of the initial 
probability, which travels within the accelerator mode. This quantity is called the survival 
probability and shown in Fig. [3j Fitting the long time decay of this probability with an 
exponential function gives us an estimate of the tunneling rate, see an example in Fig. |3j In 
our numerical simulations we take the following momentum range: 

[n(t) - 15 , n(t) + 15], n(t)=n + vt, (31) 

where no and v are given respectively by ( ]A3j) and ( 129|) . We checked that the peak traveling 
ballistically has a width less than 30 (in two photonic recoil units for the experiment, see, 



e.g., 11]) with our choice of parameters. 

We note that the right hand side in Eq. ( 130]) is defined up to a proportionality factor. 
Moreover, it crucially depends on k% and k 2 , because island sizes and tunneling rates vary 
when the kicking strength is changed. In our simulations the proportionality factor was 
chosen such as to fit the earlier regime of exponential decay (approximately between 100 
and 10 4 kicks in Fig. |3]for instance). In our numerical computations the initial state is given 
by (|52) with a width a 2 = 0.25. As the fidelity is a wildly oscillating function, it is averaged 
over 200 kicks in order to clearly expose the mean long time behavior. Results are shown in 
the Figs. H]l5] and [6] for different sets of parameters and will be discussed in the following. 
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number of kicks 

FIG. 3: (Color online) Comparison between the survival probability (black and blue dashed lines) 
and the fidelity (red full line), r = 5.86, e = t-2-k, rj = 0.01579r, /3 = 0.48984326, k 2 = 0.8vr, fe a = 
0.7-7T. The upper blue dashed line corresponds to k\ and gives the tunneling rates Ti = 5.1 x 10 -4 
whereas the lower black dashed line corresponds to k 2 and gives T 2 = 4.4 x 10~ 5 . The dot-dashed 
line corresponds to an exponential fit with a decay rate T = 3.5 x 10~ 5 . 




10000 20(100 30000 40000 50000 10000 20000 30000 40000 50000 10000 20000 30001) 40000 50000 



number of kicks number of kicks number of kicks 

a) b) c) 

FIG. 4: (Color online) Comparison between the ansatz Eq. (|30p (red dashed line) and the smoothed 
fidelity (black solid line) at r = 5.86, e = r - 2vr, rj = 0.01579r, j3 = 0.48984326, k 2 = 0.8vr and a) 
kx = 0.7vr, b) ki = 0.72vr, c) h = 0.83vr. 

The mean behavior of the fidelity for long times is well reproduced by the ansatz (130]) when 
Ti and T 2 are quite different: then the fidelity shows successively two different decay regimes, 
which are well reproduced by the pseudo-classical ansatz (130]) . see Figs. (r x ~ 10r 2 ) and 
lib (T\ ~ 5T 2 ). On the contrary, when Ti and T 2 are close to each other, one can see only 
one decay, which is still well reproduced by ( 1301) . see Fig. H] (r x ~ T 2 ) and Fig. [5] (Ti ~ 3r 2 ). 
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10000 20000 30000 40000 50000 10000 20000 30000 40000 50000 
number of kicks number of kicks 



a) b) 
FIG. 5: (Color online) Same as in Fig. Hat r = 6.6, e = r-2vr, rj = (\/5-l)/20, (3 = 0.123456789, 
&2 = 2.5 + tt] and a) k\ = 2 + rrj, b) ki = 3.5 + tt\. 




'() 10000 20000 30000 40000 50000 0001 o 10000 20000 30000 40000 50000 0001 o 10000 20000 30000 40000 50000 

number of kicks number of kicks number of kicks 



a) b) c) 

FIG. 6: (Color online) Same as Fig. U at e = -0.5, r = 4vr + e, rj = 0.001, (3 = 0.123456789, 
k 2 = (1.35 + rr])/\e\ and a) ki = (1 + rrj)/\e\, b) h = (2 + rj])/\e\, c) h = (2.2 + T V )/\e\. 

For sake of comparison we are showing the same plots for I = 2 in Fig. O The quantum 
resonance is then r = 4n. It can be seen that the agreement is not so good in the latter 
case. One reason for this may be that higher order e— classical phase space structures have 
a larger area hence may play a more important role, making estimates of the various areas 
in (1301) more problematic. Both e— classical phase spaces corresponding to k\ and hi are 
displayed in Fig. [7] corresponding to Fig. and Fig. Eb. It is clear that the shape and the 
overlap between the two islands are qualitatively different in these two situations. 
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-3-2-10 1 2 3 



FIG. 7: Phase spaces corresponding to two sets of parameters. Left column: r = 5.86, e = r — 2ir, 
T) = 0.01579t. Top left: k = 0.7vr. Bottom left: k = 0.8tt. Right column: e = -0.5, r = 4vr + e, 
r] = 0.001. Top right: k = (1.35 + rrj)/\e\. Bottom right: k = (2.2 + Trj)/\e\. 

V. CONCLUSION 



We have presented a theoretical analysis of the temporal dependence of fidelity for the 
quantum kicked rotor subject to an additional gravity field. The two major results concern 
the dynamics of this system at principal quantum resonances, i.e., at kicking periods r = 2-kI 
(I integer), and close to these resonances. In the former case, we arrive at analytical estimates 
for the decay of fidelity which is highly sensitive to the arithmetic properties of the gravity 
parameter. Close to a resonance we have used the e— semiclassical method in order to 
describe the long time asymptotics of the fidelity. The ansatz ( 130]) based on semiclassical 
densities gives a good description of the long time behavior of the quantum fidelity for both 
similar and different tunneling rates of the two compared non- dispersive wave packets 32 1 
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centered at the accelerator mode islands in phase space lg|. This result highlights once 
more the utility of e— classics in describing the quantum evolution of the kicked rotor and 
its variants. 
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Appendix A: Relevant range for the parameters. 

It was observed numerically, when plotting the e— classical phase portrait, that higher 
order nonlinear resonances play a bigger role when I is increased. For this reason we restrict 
in this paper mainly to I = 1. In the experiments, see e.g. 241 ]. k typically runs from 
0.37T ~ 0.94 to 1.57T ~ 4.71. In order to see an accelerator mode one needs a stable fixed 
point of the classical map. Following (1271) . for a given k one has a fixed point when k > tt]. 
Taking into account experimental constraints leads us to choose e such that: 

0.4<|e|<l. (Al) 

For the quantum fidelity the initial state is chosen to be a Gaussian state: 

-(n-n ) 2 /4a 2 +ine 

<*(°)> = (27ra2) i/4 » (A2) 



with 

"vrZ + r (/? + !)] , (A3) 



tt] 2jm 1 

sin 6*0 = -, n 



ek \e\ e 

where m is an arbitrary integer. This initial state is centered on an accelerator mode for 
which J = 0. This mode is travelling (in momentum space) at the speed ( )29|) . see e.g. js]. 
The fidelity is computed by applying successively the operators © for two different values 
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of the kicking strength: k\ and k-i- For each set of parameters we keep A; 2 fixed and vary 
k\. Our initial state is always chosen such as to follow the accelerator mode attached to ki\ 
once t (and e), rj and /3 are fixed, the initial state is centered in momentum space around 
n defined by ( 1A3I) . 

Appendix B: Estimating the size of the islands in the e— classical phase space 

First we choose a set of parameters r and rj for which we can have accelerator modes. 
Then we vary k within the range of existence of these modes. One way to see this range is to 
compute numerically the area of the stable island in the e— classical phase space, see Fig. |HJ 
The area follows a bell shape as a function of k. In Fig. [8] some jumps are also visible (see 

5 
4 

w: 3 

<2 
1 

°0 4 8 12 

k 

FIG. 8: (Color online) Area of the stable island in the e— classical phase space as a function 
of the kicking strength for three sets of parameters. Black full line: r = 5.86, e = r — 2ir ~ 
-0.42, 7] = 0.01579T. Red dashed line: r = 6.6, e = r - 2vr ~ 0.32, rj = 0/10. Blue dash-dotted 
line: e = -0.5, r = 4vr + e, rj = 0.001. 

e.g. for the blue dash-dotted line, between k ~ 4.424 and k ~ 4.624). We believe that this 
is due to the lack of precision when determining the island boundary and/or the breaking of 
outermost tori and their remnants (cantori). The size of the stable island in the e— classical 
phase space is computed by starting a fairly small number of trajectories outside the island. 
These are typically iterated for a long time (10 8 kicks). Then we move to polar coordinates 
((f,I) centered at the fixed point under interest (33|. The boundary of the island is then 
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determined by the curve !(</?) defined in the following way. Using a grid of thickness 5tp 
along the <p axis, the boundary Ij = I((pi) is given by: 

Jj = min {Ij, (ifj, Ij) iterated points and ipi < <pj < <fi + 5ip} (Bl) 

The size of the island is given by the area under this curve. Numerically it is computed via 
a Riemannian sum. 

The measures of the different sets in ( 1301) are computed by propagating a cloud of 10 4 
classical points. The initial points are distributed following normal distributions with mean 
and width (6o,<Je) in the 9 direction and (Jq,<jj) in the J direction. The areas in (|30|) 
reached stationary values after ca 500 kicks. The measures needed in (1301 are simply given 
by the number of points sitting on one or both of the stable islands associated to k\ and &2, 
respectively. 
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